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David L. Cheung * 1 and Friederike Schmid 1 

1 Theoretische Physik, Universitdt Bielefeld, 33615 Bielefeld, Germany 

The effect of surface roughness on the structure of liquid crystalline fluids near solid substrates 
is studied by Monte Carlo simulations. The liquid crystal is modelled as a fluid of soft ellipsoidal 
molecules and the substrate is modelled as a hard wall that excludes the centres of mass of the 
fluid molecules. Surface roughness is introduced by embedding a number of molecules with random 
positions and orientations within the wall. It is found that the density and order near the wall 
are reduced as the wall becomes rougher (i.e. the number of embedded molecules is increased). 
Anchoring coefficients are determined from fluctuations in the reciprocal space order tensor. It is 
found that the anchoring strength decreases with increasing surface roughness. 



I. INTRODUCTION 

The interaction between liquid crystalline (LC) fluids 
and solid surfaces has attracted much interest 1 . The 
presence of the surface breaks the symmetry of the LC 
phase. As well as being intrinsically interesting this is 
technologically important - many applications of liquid 
crystals depend on the interaction between the fluid and 
an external field, strongly influenced by coupling with 
external surfaces. 

Most previous studies of LC surface anchoring have 
assumed that the surface is homogenous. Two models 
are commonly used. In the first the wall is modelled 
by a perfect crystalline array 2 . The second, more coarse 
grained model, uses an external potential function that 
depends only on the distance from the wall 3,4 . While at- 
tractive from a theoretical standpoint, it has long been 
recognised that deviations from these ideal surfaces can 
affect the properties of the surface 5 . One notable ex- 
ample of this is the reduction of the order parameter 
of nematic liquid crystals at SiO surfaces^!. This con- 
trasts with measurements made on other surfaces 1 (e.g. 
rubbed polyimide) and with most simulation and theo- 
retical studies that give a higher order parameter at the 
LC-solid interface. Electron micrographs show that SiO 
surfaces are extremely rough 8 , which gives rise to the 
disordering effect of the surface. 

In this paper the structures of nematic and isotropic 
fluids near rough walls are studied. The effect of rough- 
ness is incorporated by embedding a number of molecules 
in an otherwise smooth wall. These are placed and ori- 
entated randomly. Similar models have been used for 
simple fluids 9,10 and it is hoped that this simple model 
may give insights into the behaviour of molecular fluids 
near rough or porous surfaces. Two aspects of the ef- 
fect of the surface roughness on the LC fluid are studied. 
Firstly the change in the structure of the fluid was ex- 
amined. Secondly the effect of surface roughness on the 
anchoring properties of the LC. The contribution of this 



surface anchoring to the free energy is often taken to be 
of the Rapini-Popoular formAi 



F sur f = W sin 2 (6-60) 



(1) 
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where — Oq is the angle between the director at the 
surface and the surface's 'easy-axis'. W is the surface 
anchoring coefficient. This depends on both the prop- 
erties of the bulk liquid crystal and on the interaction 
between the liquid crystal and the surface, so may be 
expected to vary with surface roughness. As this is a 
key property in applications of liquid crystals it would 
be interesting to see how this is affected by changes in 
the surface morphology. 

This paper is organised as follows. Details of the sim- 
ulation, including the method used for calculating the 
anchoring coefficient, are given in the next section. The 
structure of the fluid confined between rough walls is 
given in Sec. IIIII while results for the anchoring coeffi- 
cient are presented in Sec. II VI Finally some concluding 
remarks are given in Sec. 



II. SIMULATION 
A. Simulated Systems 

In order to simulate large systems, a simple inter molec- 
ular potential is used. This models the fluid as a system 
of soft ellipsoidal molecules interacting through a simpli- 
fied version of the popular Gay-Berne (GB) potential^. 
In particular this has two major simplifications. First 
the orientation dependence of the energy parameter is 
suppressed. Secondly the potential is cut off and shifted 
at the potential minima. These changes lead to a much 
simplier phase diagram than the GB potential, showing 
only nematic and isotropic phases, closer to the phase be- 
haviour of the hard ellipsoid 13 or hard gaussian overlap^ 
potentials. This potential is also more computationally 
efficient than the full GB potential. 

The interaction between two molecules i and j, with 
positions and and orientations and uj is given 



2 



by 



,Ui,Uj) = | 



[p- 12 - p~ 6 ] +1, p< 2 1 / 6 
, otherwise 



where eo is the energy unit, r^- = — r^, and 



(2) 



(3) 



= rij/rij, and do is the a(r^, u^, Uj) is the 



shape function given by 
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,u,) 2 



1 - xui-u^- 



-1/2 



(4) 



This approximates the contact distance between two el- 
lipsoids. In Eq. 21x = (ft 2 — l)/(ft 2 + 1) is the anisotropy 
parameter, where ft is the elongation (for the molecules 
studied here ft = 3). 

The wall is represented by a hard core potential act- 
ing upon the centres of mass of the molecules. Previous 
studies have shown that this gives rise to homeotropic 
alignment at the wall 16 . Roughness is introduced by em- 
bedding a number of molecules, N w in the wall. These 
were given random positions and orientations which were 
kept fixed during the simulations. While generating these 
surface configurations interactions between the surface 
molecules were ignored, thus these molecules may over- 
lap. It should be noted that these molecules do not cor- 
respond to real molecules, rather they are used as a con- 
vient way of introducing inhomogenity into the wall. The 
roughness of the wall was characterised by the surface 
density of these embedded molecules E = N w /A. Some 
example wall structures are shown in Fig. To ensure 
some sampling of surface configurations three different 
surfaces were studied for each pair of p and E. 

Simulations were performed at two average densities, 
p = 0.314 and p = 0.30. For the higher density the 
fluid confined between smooth walls was nematic, while 
it is isotropic for the lower density. The simulated sys- 
tems were composed of 1200 fluid molecules and up to 
63 molecules embedded in each wall. Throughout this 
reduced units defined by the molecular width cjq and the 
energy unit eo are used. A reduced temperature of 0.5 
was used for both densities. 



B. Simulation observables 




FIG. 1: (Color online) Example rough wall configurations for 
(a) E 0.2 and (b) E 0.4. 



where is the orientation of the zth molecule and 8 a p 
is the Kronecker delta function. It may also be infor- 
mative to consider the order parameter in the cell bulk 
and near the surface, S bulk and S sur f . S bulk is calcu- 
lated for molecules within the region l z /4 < z < 3Z 2 /4, 
while S sur f is calculated for molecules within 1 <jo of the 
surface. 

The distribution of molecules in the simulation cell can 
be described by the density profile p(z). To describe 
the ordering through the cell, the ordering tensor Eq. 
can be calculated throughout the cell. Diagonalising 
this gives the order parameter profiles (q+(z), qo(z), and 
q-(z)). These can be expressed as S(z), S(z) + \S xy (z), 
and S(z) — ^S xy (z), where S(z) is the nematic order 
parameter and S xy (z) is the biaxiality parameter. 

The nematic director n(z) can be identified with the 
eigenvector of the ordering tensor corresponding to the 
largest eigenvalue. 

While the presence of layers may be deduced from the 
density profiles it may be useful to quantify the degree 
of translational order. The smectic order parameter may 
be introduced for this purpose 17,18 . This is given by 



The orientational order may be characterised by the 
usual nematic order parameter. This is given by the 
largest eigenvalue of the ordering tensor, defined as 



Pi 



N 



1 ^ / 2iriZj 

ivS exp T" 



(6) 
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a,[3 = x,y,z (5) 



where d is the layer periodicity. This is initially unknown 
and is take to be the value that maximises pi— . 



3 



C. Director fluctuations and surface anchoring 



III. FLUID STRUCTURE 



The surface anchoring coefficient is determined by the 
director fluctuation methocU2*2&. This method relates 
thermal director fluctuations in a confined geometry to 
the zenithal anchoring coefficient, in a similar manner as 
the fluctuations in a bulk LC can be related to the bulk 
elastic constantsSiiSi. The theory for this has been exten- 
sively developed elsewhere and this section will contain 
only the briefest of outlines. 

As for the bulk elastic constants the zenithal anchor- 
ing coefficient may be determined by fitting elastic the- 
ory predictions of fluctuations in the ordering tensor to 
those determined from simulations. The reciprocal space 
ordering tensor is given by 



V 



Qa/3(k) = J^Yl Git* eXP ( ik - r ^ * 



Fluctuations can be calculated from simulation 



(7) 



9\ V 
\Qap(k z )\ 2 ) = ^ 



J2®ip cos(k z Zj) 



, (8) 



The corresponding elastic theory predicts that there fluc- 
tuations are given by 19 



x 2 +e 



£\ e i("~x) - 1 



(9) 



where K33 is the bend elastic constant. q z is a wave 
vector with a discrete spectrumi2, x = QzL e , and k = 
k z L e . £ is the anchoring strength parameter 



33 



Le 
A 



(10) 



where W is the zenithal anchoring coefficient and A is the 
extrapolation length. L e is the cell thickness appearing 
in the elastic theory; this is not necessarily equal to the 
simulation cell thickness, L z . In fitting the elastic theory 
to simulation profiles L e and £ are the fitting parameters. 
K33 has been determined from simulation for a nearby 
state point 23 (p = 0.30). While this value (if 33 = 1-48) 
is likely to be too large for some of the systems studied 
here, this should be sufficient for a qualitative study. 



A. High Density Fluid 

The density profiles for the high density fluid are shown 
in Fig. Et 1 The effect of the wall roughness is most ap- 
parent near the wall. Here the density near the surface 
decreases with increasing E. This is caused by the de- 
crease in available volume near the wall due to the em- 
bedded molecules. Values of the density near the wall 
are presented in Tab. [I] The surface density falls from 
0.72 for the smooth wall to 0.34 for the rough wall with 
E = 0.4. 

Another noticeable difference is that the second peak 
(at z = 2.8 for the plain wall) becomes broader. This 
arises from the surface disorder disturbing the layer struc- 
ture and has been observed in simulations of Lennard- 
Jones fluidslfi-. This can more clearly be seen in the in- 
set, which shows the detail of the density profiles around 
the minima. The disruption of the translational order- 
ing caused by the embedded molecules can be seen by 
considering the smectic order parameter (Eq. EJ). Values 
for these are presented in Tab. As can be seen p\ 
markedly decreases with increasing grafting density, as 
would be expected for increasing translational disorder. 

Far from the wall the profiles all tend to a constant 
values, indicating a layer of bulk fluid. The density of 
this layer increases slightly with increasing grafting den- 
sity. This arises as the embedded molecules exclude fluid 
molecules from regions near the wall, increasing the num- 
ber of molecules in the cell bulk. This is a consequence 
of having a fixed cell size and may be avoided by using 
NpT simulations. Quantitively this can be seen by ex- 
amining the densities in the cell bulk. Values for this are 
presented in Tab. Q The density in the bulk of the cell 
goes from 0.29 for the smooth wall to 0.31 for the highest 
grafting densities. 

Figure Efb) shows the order parameter profiles for dif- 
ferent values of E. As can be seen the value of the or- 
der parameter at the wall is lower for higher grafting 
densities. This is caused by the disorientating effect of 
the embedded molecules. This disorientating effect also 
leads to a deeper minima. For E > 0.2 this leads to 
a small layer (approximately 1 molecular width thick) 
of almost isotropic fluid. The position of this minima 
moves closer to the wall with increasing surface rough- 
ness. For the smooth wall this minima is at approxi- 
mately z = 2, while for the highest grafting densities it 
appears at about z = 1.1. Again this is attributable to 
the disruption in the surface induced layering. As for 
p(z) the second peak becomes broader with increasing 
E. Finally, as can be seen from Tab. Q]the bulk order 
parameter S bulk increases with increasing E. This is a 
consequence of the increasing density in the centre of the 
cell due to the excluded volume effect of the embedded 
molecules. It is noticeable that for E > 0.2 S bulk becomes 
larger than S surf . 

The biaxiality profiles are shown in Fig. Etc). For the 
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FIG. 2: (a) Density profiles for the high density fluid near 
rough walls. The density profile for grafting density S = 
is shown by the solid line, for E = 0.1 dotted line, E = 
0.2 dashed line, E = 0.3 long dashed line, and E = 0.4 the 
dashed dotted line. Inset shows the density profiles around 
the minima. Symbols as in main figure. 

(b) Order parameter profiles for high density fluid near rough 
walls. Symbols as in (a). 

(c) Biaxiality (S xy ) profiles for high density fluid near rough 
wall. Symbols as in (a). 

(d) z component of the director for the high density fluid. 
Symbols as in (a). 3 n z (z) profiles are shown for E = 0.4. 



smooth wall the this is essentially zero (the largest value 
is 0.04) reflecting the cylindrical symmetry around the 
z axis. However, for the rough walls there is are sizable 
peaks in the biaxiality profiles. These are stronger for 
larger values of E and are in the region of 0.5 < z < 
1.3, corresponding to regions of low order. This surface 
induced biaxiality has been seen for simulations of LCs 
near grooved surfaces 24 . 



Table 1. 

Densities and order parameters for the simulated 
systems, pbulk and p SU rf are the bulk and surface 
densities, S, S bulk and S surf are the total, bulk and 
surface order parameters, and p\ is the smectic order 
parameter. Errors in the last decimal place are in 
parenthesises. 
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E 


Psurf 


Pbulk 


b 




obulk 


Pi 


0.314 





0.72(1) 


0.286(3) 


0.60(3) 


0.84(1) 


0.53(6) 


0.14(2) 


0.314 


0.1 


0.61(1) 


0.294(2) 


0.64(2) 


0.76(2) 


0.65(3) 


0.12(2) 


0.314 


0.2 


0.48(3) 


0.304(3) 


0.67(2) 


0.65(3) 


0.72(2) 


0.10(2) 


0.314 


0.3 


0.43(5) 


0.308(4) 


0.69(2) 


0.65(9) 


0.75(2) 


0.09(2) 


0.314 


0.4 


0.33(2) 


0.307(9) 


0.66(8) 


0.58(7) 


0.73(7) 


0.07(2) 


0.300 





0.70(1) 


0.273(2) 


0.28(3) 


0.81(1) 


0.10(4) 


0.14(2) 


0.300 


0.1 


0.59(1) 


0.281(2) 


0.34(7) 


0.72(3) 


0.27(9) 


0.12(2) 


0.300 


0.2 


0.46(2) 


0.291(3) 


0.52(6) 


0.61(3) 


0.59(5) 


0.10(2) 


0.300 


0.3 


0.41(2) 


0.293(3) 


0.51(4) 


0.53(9) 


0.60(4) 


0.09(2) 


0.300 


0.4 


0.38(1) 


0.300(3) 


0.59(6) 


0.54(5) 


0.69(3) 


0.07(2) 



Figure EJd) shows the z component of the director for 
each E. In the cell bulk the director is essentially parallel 
to the z axis. For E > 0.2 there is a tilt away from the 
z axis at about the position of the order minima. As 
may be expected this is most pronounced for the E = 0.4 
wall. In Fig. [2]i the profiles for each of the E = 0.4 walls 
are shown separately. It can be seen that the size of this 
tilt differs strongly for different wall configurations (for 
the largest the tile angle is approximately 79°). For the 
larger tilt angles this propagates into the bulk of the fluid 
leading to a director tilted up to 16° from the z-axis. It 
is not clear how a randomly generated wall gives rise to 
a titled configuration in the bulk. Similar behaviour has 
been seen in a recent study of a LC near a planar wall 
with perpendicularly grafted rods» In that case the 
bulk tilt was caused by the competition between the wall 
(which promoted planar alignment) and the embedded 
molecules. As it appears only for a subset of the walls 
studied here it would be desirable to consider further wall 
configurations. 

The previous discussion may be illuminated by exami- 
nation of simulation configurations. Figure |3| shows con- 
figurations of for E = 0.0 (smooth wall), E = 0.2, and 
E = 0.4. The disordering effect of the rough wall can be 
seen in the first and second layers (left and right most 
pictures). However, the fluid between these two layers 
shows the most noticeable change with increasing E. For 
the smooth wall the molecules in this region are still well 
ordered parallel to the z-axis. With increasing E the 
molecules in this region become increasingly disordered. 
This gives rise to the deeper minima seen in the order pa- 
rameter profile (Fig. Efb)). Additionally it can be seen 
that many of the molecules lie in the xy plane, giving rise 
to the biaxiality peak and the tilt of the director away 
from the z axis. This behaviour is similar to that seen in 
simulations of smectic liquid crystals 26 where molecules 
in the region between the layers are seen to align either 
parallel or normal to the layers. These planar oriented 
molecules possibly give rise to the bulk tilt seen in some 
cases. Finally the number of molecules in this region 
visibly increases with E. 
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FIG. 3: (Color online) Simulation configurations showing 
molecules within 2.5 ao of the surface for (a) E = 0, (b) 
£ = 0.2, and (c) E = 0.4. For each E the left most picture 
shows molecules with < z < 0.5, the centre pictures shows 
0.5 < z < 1.5, and the right most shows 1.5 < z < 2.5. 




FIG. 4: (a) Density profiles for the low density fluid near 
rough walls. The density profile for grafting density E = 
is shown by the solid line, for E = 0.1 dotted line, E = 0.2 
dashed line, E = 0.3 long dashed line, and E = 0.4 the dashed 
dotted line. 

(b) Order parameter profiles for the low density fluid. Sym- 
bols as in (a). 
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B. Low Density Fluid 

Here the density and order parameter profiles for the 
low density system are discussed. For the smooth wall the 
density in the bulk of the cell is 0.27 (Tab. [Q, just below 
the isotropic-nematic transition density for this system 
(pi-N — 0.287). As the density in the cell bulk increases 
with E, for E > 0.2 the fluid in the cell bulk is nematic 
rather than isotropic. 

The density profiles for the low density fluid are shown 
in Fig. Ufa). The changes in the density profile with in- 
creasing E are similar to those in the high density system 
- the density at contact decreases with E and the second 
peak becomes more diffuse. Again this can be gleaned 
from the decrease in the value of the smectic order pa- 
rameter with E (Tab. ^J. It is interesting to note that 
the values of p\ obtained in this system are very simi- 
lar to those for the higher density system, indicating the 
similarity in the structure of both systems. 

Shown in Fig. Efb) are the order parameter profiles. 
As in the high density fluid the value of the order pa- 
rameter at the wall decreases as E increases. The order 
parameter profile also shows a deeper minima with in- 
creasing surface roughness. It is noticeable that even in 
this lower density case there is not an appreciable layer 
of isotropic fluid between the wall and bulk fluid. This 
has been predicted to happen near rough walls as a con- 
sequence of the competition between the bulk director 



FIG. 5: Order tensor fluctuations (normalised by cell volume) 
as a function of wavevector for (a) high density and (b) low 
density fluids. In both graphs the simulation data is denoted 
by symbols (circles E = 0, squares E = 0.1, diamonds E = 
0.2, triangles E = 0.3, and crosses E = 0.4) and the elastic 
theory data is shown by lines (continuous line E = 0, dotted 
line E = 0.1, dashed line E = 0.2, long dashed line E = 0.3, 
and dot dashed line E = 0.4). The order tensor fluctuations 
for E = 0.0 are shown only in (a). 

and the local boundary conditions 27 . 

IV. SURFACE ANCHORING 

Shown in Fig. [51 are the order tensor fluctuations as a 
function of wavevector. As can be there is good agree- 
ment between the simulation and elastic theory curves, 
especially for small k z . 

The fitted values for the anchoring coefficients are 
given in Tab. along with values of the extrapolation 
length A and the surface anchoring coefficient W. For 
both bulk densities ( tends to decrease with increasing 
E. 

Table 2. 

Fitting data for the order tensor fluctuations (Fig. 0. £ 
is the anchoring coefficient, and L e is the elastic theory 
cell width, which appear in Eq. A = L e /( is the 
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extrapolation length and W = if 33 /A is the surface 
anchoring strength. 



p 


S 


c 




A 


W 


0.314 


0.0 


5.62 


16.00 


2.85 


0.52 


0.314 


0.1 


6.51 


28.39 


4.36 


0.34 


0.314 


0.2 


5.64 


29.36 


5.21 


0.28 


0.314 


0.3 


4.94 


26.80 


5.43 


0.27 


0.314 


0.4 


4.55 


27.78 


6.11 


0.24 


0.30 


0.1 


3.05 


18.32 


6.01 


0.24 


0.30 


0.2 


2.88 


24.60 


8.54 


0.17 


0.30 


0.3 


2.58 


23.85 


9.24 


0.16 


0.30 


0.4 


2.70 


25.44 


9.42 


0.16 



The behaviour of the elastic theory cell width, L e , with 
E deserves comment. For the high density fluid, L e for 
the smooth wall is 16.00 <jo, a few molecular lengths 
smaller than the physical cell width (L z = 24.66). This 
is similar to behaviour seen in previous simulations 19,20 
and is due to the formation of highly ordered layers in 
the vicinity of the surface. For the rough walls however, 
L e is larger than L z . This increase may be attributable 
to the rough surface breaking up the highly ordered sur- 
face layer. Thus instead of the elastic theory boundary 
conditions being applied at this layer, they are applied 
closer to the wall, leading to an increase in L e . 

For both bulk densities the extrapolation length in- 
creases and the anchoring coefficient W decreases with 
E. Thus, as may be intuitively expected, anchoring on 
rough surfaces is weaker than on smooth surfaces. 

V. SUMMARY 



corresponds to an increasing surface roughness. The sim- 
ulations were performed at two bulk densities. For the 
higher density the fluid in the bulk of the cell is nematic 
for the smooth wall case, while for the lower density it is 
isotropic. 

At both densities studied the effect of increasing sur- 
face roughness is similar. Both the density and order 
parameter in the region near the wall decrease as the 
number of embedded molecules increases. The decrease 
in the density arises from the excluded volume of the 
embedded molecules, while the decrease in the order can 
be attributed to the disorientating effect of the randomly 
orientated molecules in the wall. The rough walls also 
act to smear out the secondary peaks in the density and 
order parameter profiles as the embedded molecules give 
anchoring points at positions other than at the wall sur- 
face. One side effect of the wall roughness is an increase 
in the density and order parameter in the centre of the 
cell. 

Also studied was the effect of surface roughness on the 
surface anchoring strength. For both systems the an- 
choring was found to become progressively weaker with 
increasing surface roughness. 

A number of possible avenues for future work are pos- 
sible. Calculation of the anchoring coefficient via alterna- 
tive methodsi£i2£ would be useful. As the formation of a 
highly ordered layer at the surface is commonly held to be 
important for the growth of order in confined liquid crys- 
tals, it may be interesting to investigate the effect of wall 
roughness on the phase behaviour of the confined fluid 17 . 
Integral equation^ 9 - or density function theories^ have 
been applied to similar systems of simple fluids and ap- 
propriate generalisations to molecular fluids should also 
prove useful. 



In this paper results of Monte Carlo simulations for 
a confined fluid of ellipsoidal molecules have been per- 
formed. The effect of surface roughness on the structure 
of the fluid has been examined. Roughness was intro- 
duced by embedding a number of molecules, with random 
positions and orientations, in otherwise smooth walls. In- 
creasing the number of molecules embedded in the wall 
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